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I  INTRODUCTION 


A.  Background  ^  ' 

In  the  last  ten  years,  high-speed  (10-  to  100-GHz)  devices  have  been  tested  in  myriad 
operating  modes  and  geometries.  The  goal  of  these  efforts  have  been  to  develop  efficient 
discrete  elements  that  can  be  readily  fabricated  into  reliable  systems  with  high  production  yield. 
While  the  tried  and  true  silicon  semiconductor  continues  to  be  the  workhorse  of  the  industry, 
its  relatively  high  effective  mass  begins,  makes  it  less  attractive  that  some  of  the  III-V  com¬ 
pounds  (e.g.,  GaAs,  InP,  and  certain  alloys)  once  design  frequencies  are  high. 


Moreover,  many  of  the  new  device  designs  call  for  heterojunctions  (e.g.,  tunneling  junc¬ 
tions  in  hot  electron  transistors),  high  mobility  channels  (e.g.,  quantum  well  structures)  and 
other  special  characteristics.  The  question  of  optimum  materials  selection  for  these  multima¬ 
terial  structures  is  nontrivial.'  The  advent  of  such  growth  methods  as  molecular  beam  epitaxy 
(MBE)  and  metal  organic  chemical  vapor  deposition  (MOCVD)  have  made  it  possible  to 
prepare  multilayer  thin  films  of  different  materials  for  each  different  part  of  a  device  (e.g.,  the 
emitter,  collector,  and  base  of  a  transistor).  There  are  general  guidelines  within  which  materi¬ 
als  are  compatible.  Most  of  the  semiconductors  of  interest  fall  into  the  diamond  or  zincblende 
crystal  structure  groups,  which  can  be  grown  compatibly  on  dissimilar  substrates  without 
device-degrading  defects  provided  that  the  lattice  constants  of  the  substrate  and  epitaxial  layers 
are  not  too  different. 


-  In  this  report,  we  address  only  a  subset  of  questions  bearing  on  materials  selection;  how¬ 
ever,  these  questions  arc  central  to  high-speed  performance.  The  major  question  is,  which 
materials  offer  the  best  set  of  performance-limiting  parameters  for  the  device’s  active  transport 
region  (usually  the  base)?  The  number  of  possible  materials  from  which  to  select  is  enormous: 
the  compounds  and  alloys  made  from  the  cations,  aluminum,  gallium,  indium,  and  the  anions 
phosphorous,  arsenic,  and  antimony.  Thus,  predictions  of  an  accurate  theory  can  serve  to 
reduce  greatly  the  number  of  experiments  needed  to  optimize  performance. 


B.  Device  Operating  Modes 

The  answer  to  our  question  differs  for  devices  intended  to  function  in  grossly  different 
operating  modes.  The  modes  can  be  separated  into  three  general  classes;  drift,  ballistic,  and 
tunneling.  The  distinction  between  the  classes  is  determined  by  the  length  of  the  active  region 
of  the  device. 

•  If  the  device  is  long  compared  to  the  majority  carrier  velocity  mean  free  path,  then  it 
functions  in  the  drift  mode,  in  which  carriers  make  many  collisions  as  they  transport 
across  the  device.  The  mean  free  path  reduces  to  the  momentum  mean  free  path  only  if 
the  effective  mass  is  a  constant  for  all  the  populated  states.  Generally,  the  lower  size 
limit  for  the  drift  mode  is  a  few  tenths  of  a  micrometer. 
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.  Ballistic  mode  devices  are  those  whose  scale  size  falls  between  the  velocity  mean  free 
path,  and  the  tunneling  thickness,  generally  tens  of  angstroms. 

.  The  tunneling  mode  devices  are  the  smallest  ones. 

Ballistic  mode  devices  have  most  of  the  carriers  transport  through  the  active  region  with  few  or 
no  collisions.  Most  of  this  report  is  devoted  to  materials  choices  for  ballistic  mode  devices, 
and  nothing  more  will  be  said  about  the  tunneling  modes. 

The  differing  criteria  for  materials  choices  in  the  drift  and  ballistic  modes  lead  to  distinct 
optimum  ch'-'ces. 

.  All  of  the  drift  mode  device  performance-limiting-parameters  can  be  expressed  in  terms 
of  features  of  the  drift  velocity-field  curves:  At  low  fields,  these  curves  increase  linearly 
with  field  (the  slope  is  the  low  field  mobility).  For  some  materials  (e.g.,  silicon  and 
other  indirect  bandgap  semiconductors),  these  curves  then  saturate;  in  others  (GaAs, 
InP,  and  most  direct-bandgap  semiconductors),  they  first  reach  a  peak,  then  saturate 
after  descending  through  a  valley. 

.  The  ballistic  mode  places  a  premium  on  the  speed  with  which  the  carrier  moves,  and 
the  fraction  of  electrons  reaching  the  collector  without  colliding.  The  speed  and  device 
dimensions  determine  the  upper  cutoff  frequency  and  the  ballistic  fraction  is  related  to 
efficiency. 


C.  Significance  of  Band  Structure 

A  reliable  prediction  of  the  device-performance-limiting  parameters  of  high-electron- 
dominated  high-speed  structures  requires  the  use  of  accurate  band  structures,  rather  than 
effective-mass  approximations.  This  is  because  the  carriers  occupy  highly  excited  states  in 
which  the  curvature  of  the  bands  ek  changes  rapidly  with  crystal  momentum  k.  This  affects 
not  only  the  effective  mass,  but  also  the  group  velocity  and  the  density  of  states. 

The  density  of  states  is  particularly  important  in  calculating  scattering  rate.  It  has  been 
traditional  to  calculate  scattering  rates  from  effective-mass  approximations.  Coupling  constants 
for  some  scattering  mechanisms  (e.g.,  electron  to  polar  optical  phonons)  are  reliably  calculated, 
while  others  (notably  electron  to  deformation-potential-coupled  longitudinal  acoustic  phonons) 
are  not.  [We  here  restrict  our  discussions  to  electrons,  although  some  light-hole  band  devices 
have  been  suggested.]  Hence,  it  is  common  practice  to  use  effective-mass  band  structures  for 
the  central  T-vallcy  and  satellite  L-valley  when  calculating  the  polar  phonon  impurity  and 
electron-electron  contributions  to  intervalley  scattering,  then  to  choose  the  deformation  potential 
coupling  constant  such  that  the  net  scattering  rate  fits  an  experiment  (e.g.,  four-wave  optical 
mixing).  The  use  of  effective-mass  band  structures,  which  produce  densities  of  final  states  that 
are  too  small,  results  in  underestimates  of  polar  phonon  and  the  impurity  scattering.  As  a 
consequence,  the  deformation  potential  coupling  constant  needed  to  fit  experimental  data  is 
larger  than  its  true  value.  This  procedure  distorts  the  relative  importance  of  the  deformation 
potential  and  other  mechanisms.  When  the  various  mechanisms  are  applied  to  situations 
differing  from  those  of  the  experimental  data  to  which  they  were  fit,  this  procedure  can  intro¬ 
duce  serious  errors.  Thus,  accurate  band  structures  are  required  to  make  reliable  predictions 
about  hot  electron  transport  in  any  of  the  operating  modes. 
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It  is  also  desirable  to  calculate  all  the  coupling  constants  so  as  to  eliminate  the  need  for 
direct-scattering-related  experimental  information.  The  electron-scattering  mechanisms  we  have 
treated  this  way  are  longitudinal  and  transverse  for  both  optical  and  acoustic  phonons,  screened 
charged  impurity,  plasma,  and  alloy.  Alloy  band  structures  and  scattering  rates  are  calculated 
together  in  the  coherent  potential  approximation  (CPA)  formalism. 


D.  Organization  of  the  Report 

Section  II  summarizes  calculation  of  scattering  rates  resulting  from  various  scattering 
mechanisms.  We  discuss  our  choice  of  materials  for  ballistic  transport  in  Section  III  and  that 
for  drift  transport  in  Section  IV.  Because  most  of  the  results  have  appeared  in  the  form  of 
journal  publication,  the  details  are  not  given.  Conclusions  are  given  in  Section  V.  A  Bibliog¬ 
raphy  of  journal  publications  and  conference  papers  resulting  from  this  36-month  research  fol¬ 
lows  the  list  of  references. 


n  CALCULATION  OF  SCATTERING  RATES 
FROM  VARIOUS  MECHANISMS 


The  full  band  structure  is  used  in  the  calculation  of  scattering  rates  resulting  from  ionized 
impurity,  longitudinal  optic  (LO)  phonons,  deformation-potential-coupled  intervalley  scattering, 
electron-electron  scattering,  and  (in  the  case  of  alloys)  disorder  scattering.  The  inclusion  of 
full  band  structure  makes  the  Brillouin  zone  (BZ)  sum  more  time-consuming,  because  the  sum 
must  be  carried  out  over  the  entire  BZ,  instead  of  the  usual  l/48th  of  BZ.  All  coupling  con¬ 
stants  necessary  to  evaluate  the  scattering  rates  are  obtained  from  the  underlying  band  struc¬ 
ture.  Because  some  of  these  results  have  been  described  in  the  quarterly  reports,  the  procedure 
is  discussed  here  in  brief.  Details  can  be  obtained  from  our  published  papers. 

Our  tight-binding  (TB)  Hamiltonian  H0(k)  includes  all  long-range  interactions  (Krish- 
namurthy,  Sher,  and  Chen,  1986).  Our  basis  contains  the  Gaussian  orbitals  of  the  type  lodj> 
(a  can  be  s,  pv  py,  or  p,)  for  each  sublattice,  j,  and  each  cell,  1.  The  Bloch  orbitals  I  odcj>  are 
expanded  in  terms  of  large  number  of  plane  waves  I  It  +  <j  >  (d  is  reciprocal  lattice  vector) 
and  the  H0(ic)  is  obtained  by  using  pseudopotentials.  The  band  structure  is  fine-tuned  to  exper¬ 
imental  values  by  adding  an  adjustable  Hamiltonian  Hj(1c).  The  expansion  coefficients  A(IJ, 
which  are  eigenvectors  of  H0(lc)  +  H,(k5,  are  used  to  obtain  the  wave  function  required  in  the 
calculation  of  various  scattering  rate. 

A.  Ionized  Impurities 

The  potential  for  ionized  impurities  scattering  is  given  by: 

*  ' 

vei=  + 

1>V  J  1 

where 

X  =  [(4tc  n,  e2)/(tcJcBT)]1/2  , 

and  n[  is  the  impurity  density,  N,/V,  tc*  is  the  infinite  frequency  dielectric  constant,  kB  is 
Boltzmann’s  constant,  e  is  the  electron  charge,  and  T  is  the  lattice  temperature  in  Kelvins.  For 
scattering  from  k'  to  k, 
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where  a  is  the  lattice  constant  in  angstroms,  n;  is  in  1020  m  3,  and  ~h  is  Planck’s  constant 
divided  by  2tr.  In  these  expressions,  k  and  are  in  units  of  (27t/a)  and  eV,  respectively. 


B.  Phonons  in  Compounds 

The  potential  for  scattering  from  LO  phonons  in  compounds  is  given  by: 


where  aq  ,  aq+  arc  phonon  annihilation  and  creation  operators,  ct>i  is  the  longitudinal  phonon  fre¬ 
quency,  V  is  the  crystal  volume,  and  Kq  is  the  low-frequency  dielectric  constant.  For  a  scatter¬ 
ing  from  k'  to  k, 

c  Nq0fph  Pi(Xic')  (absorption) 

5Dh(K,K  )  -  . 
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In  the  calculation  of  Pi,  P2,  and  fph  terms,  k  and  are  in  (2rc/a)  and  eV,  respectively. 


C.  Intervalley  Scattering — Deformation  Coupled 

While  the  BZ  sum  of  I1Ic,kv)»  PiOcjc'),  P20c,ic')  includes  that  portion  of  intervalley  scatter 
ing  mediated  by  LO  phonons  and  impurities,  there  is  additional  intervalley  scattering  arising 
from  the  displacement  of  atoms  from  their  regular  sites.  This  deformation-potential-coupled 
intervalley  scattering  is  believed  to  be  very  important  at  high  electric  fields.  However,  the 
strength  of  this  scattering  depends  sensitively  on  the  less-known  deformation  potential  constant 

D.  We  had  made  an  ab  initio  calculation  of  this  D  from  the  underlying  band  structure  (Krish- 
namurthy,  Sher,  and  Chen,  1988c).  Calculated  scattering  rates  successfully  explain  characteris¬ 
tics  of  GaAs  HET.  The  scattering  rate  from  1?  in  T-valley  to  1c  in  L  valley,  for  example,  is 
given  by: 


siv(ic,ic') 


flv  P3(tic0 
flv  P40tf ') 


(absotption) 

(emission) 


P3(kV)  = 


P4(£k')  = 


'-'k 


1  1  k 

f,v  =  1.332  x  105 


D2 


a 3  p  <7tco0) 


x  1012  s  1 


All  1c  and  are  in  units  of  (2rc/a)  and  eV,  respectively.  D  is  in  108  eV/cm,  a  is  in  angstroms, 
p  is  in  g/cm3,  and  7jco0’  the  energy  phonons  necessary  for  intervalley  transfer,  is  in  meV. 


D.  Alloy  Disorder  Scattering 

The  electron  transport  in  alloys,  calculated  in  CPA,  is  affected  by  disorder  in  two  ways. 
The  real  part  of  self-energy  changes  the  relationship  of  e  to  k  end  thus  changes  the  density  of 
states.  The  imaginary  part  of  self-energy  introduces  a  finite  lifetime  to  these  states  and  thus 
additional  scattering  is  experienced.  In  AXB(_XC  alloys,  for  example,  the  electron  sees  different 
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potential,  depending  on  whether  an  A  or  B  atom  is  at  that  site.  If  we  denote  this  potential 
change  in  site  from  VCA  average  in  a  matrix  8V,  then  the  alloy  disorder-disorder  scattering 
potential  is: 


AV  =  —  V  <n  lSVIn>  , 

N  t 

where  n  is  the  site  index  and  N  is  total  number  of  sites.  Then,  for  a  scattering  from  t'  to  1c, 
we  have: 

SA(jc,i?)  =  i  <5cl  AV \t'>  1 2  5(&?  -  to)  , 

7i 

=  fA  h&)  . 

I2(i?,P)=]i;A;jAoi.AV»„.]2-|^r  . 

,aa'  1 

I  u  I 

fA  =  9.48  x  103  x(l  -  x)  x  1012  s_1  , 
where  k  and  Eg  are  in  units  of  (2tt/a)  and  eV,  respectively. 

E.  Phonons  in  Alloys 

In  alloys  such  as  AXB,_XC,  there  will  be  two  kinds  of  phonons,  namely  AC-type  and  BC- 
typc.  The  coupling  constant  for  interaction  with  these  phonons  will  be  strongly  concentration- 
dependent.  Moreover,  the  number  of  phonons  in  each  kind  will  also  be  concentration- 
dependent.  Hence,  polar  phonon  scattering  in  alloys  can  be  very  different  from  that  in  consti¬ 
tuent  compounds. 

The  q-dependent  dielectric  function  for  phonon-plasmon  coupled  system  is: 

(  2  2  1 
CD  -  CDi 

e(q,co)  =  — - -  i-  x(q,co)  , 

iy  -  y  J 

where  x(q,co)  in  Lindhard  form  is 
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x(q,co)  =  Xi(q.o)  +  i  x2(q.“) 


kpr 

Xi(q.co)  =  — 

q 


±  + A . 

• 

f  ,  1 

co  +  q^m 

2 

-  1 

2  4q 

qvF 

. 

. 

v  J 

. 

(0  —  qvp  +  q2/2m  1 
I  co  +  qvp  +  q2/2m  | 


co  -  q2/2m 


qvF 


-  1 


In 


|  co  -  qvF  -  q2/2m  J 
I  co  +  qvF  -  q2/2m  | 


When  q  <  2kF, 


X:(q.«)  = 


JL  ^FT  w 
2  q2  qvF 

K_  ^FT  kp 

4  q2  q 
0 


1  - 


co  -  qz/2m 
qvF 


12' 


0  <  co  <  qvF  -  qz/2m 


-q2/2m  <  co  -  qvF  <  q2/2m 
co  >  q2/2m  +  qvF 


When  q  >  2kp, 


X2(q.<o)  = 


_7t  KFT  kF 
2  q2  q 


1  - 


co  -  q2/2m 
qvF 


q2/2m  -  qvF  <  co  <  q2/2m  +  qvF 
otherwise 


with 


kpr  = 


(Thomas-Fermi  screening) 


Vep  =  i 


4kc2H 

1/2 

v  1 

A 

K.V 

.®ii 

1  + 


A 22 


M2  J 


1/2 


(a,q  e'^-  a,t  e'^ 


“12  L 


1  - 


a22 


(a2q  elq  r  -  a2q  e~,qr) 
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S(£ic')  =  1012  x  86.94 


F.  Electron-Electron  Scattering 

It  is  believed  that  in  a  heavily  doped  semiconductor  electron-electron  scattering  represents 
the  fundamental  limitation  on  the  electron  mean  free  path.  To  quantify  this  effect,  we  have 
calculated  the  velocity  relaxation  time  in  effective  mass  approximation  as  follows.  For  the 
scattering  k*  to  k', 

Sdc,icO  =  -8rc  e2  f  dca  1m  — - — -  8(&?  -  E£±  ±  H  to) 

J  £(q,d)>  ^ 

where  k'  is  k*  t  ^  and  the  ±  denotes  the  absorption  or  emission  process. 


HI  MATERIALS  CHOICES  FOR  BALLISTIC  MODE  DEVICES: 
HOT  ELECTRON  TRANSISTORS 


Both  GaAs  and  InAs,  two  important  materials,  are  direct-gap  materials;  the  central-valley 
r-minimum  is  the  bottom  of  the  conduction  band,  the  L-valley  minimum  in  the  (111)  direction 
is  next  lowest,  and  the  X-valley  in  the  (100)  direction  is  higher  still.  The  group  velocity  Vg(k) 
at  any  wave  vector  k  is  given  in  terms  of  the  band  energies  e(k)  by 

vg(i?>  =  TVe(^  ■ 

8  ti 

We  have  shown  (Krishnamurthy,  Sher,  and  Chen,  1988a,  1988b)  that,  for  the  direct-gap 
materials,  the  group  velocity  reaches  its  highest  values  for  electrons  injected  in  the  (100)  direc¬ 
tion.  Thus,  the  remaining  discussion  is  confined  to  injection  in  the  (100)  direction.  Figure  1 
shows  the  group  velocity  as  a  function  of  injected  energy  for  several  semiconductors.  The 
peak  in  vg  for  GaAs  is  vgp  =  1.2  x  108  cm/s  and  at  E  =  0.8  eV;  however,  the  curve  is  quite  flat 
over  an  extended  range  of  energies.  Only  InAs  and  the  Ga]_xInxAs  alloys  have  a  vgp  larger 
than  that  of  GaAs.  However,  because  the  band  gap  of  InAs  is  small,  electrons  injected  at  ener¬ 
gies  near  its  vgp  can  cause  a  breakdown.  The  bandgap  of  the  Gao.sInosAs  alloy  is  large 
enough  to  preclude  this  problem. 

Qearly  the  optimum  strategy  to  maximize  device  speed  is  to  inject  the  electrons  from  the 
source  at  ’heir  peak  group  velocity,  Vgp.  If  this  is  done,  any  attempt  to  accelerate  them  with 
applied  fields  is  useless,  because  any  change  in  k  will  decrease  their  velocity.  There  is  another 
concern,  however.  If  the  scattering  rate  at  vgp  is  high  enough  that  the  velocity  mean  free  path 
is  low,  and  if  few  electrons  reach  the  drain,  the  device  efficiency  will  be  low.  Thus,  a 
compromise  between  speed  and  efficiency  may  be  required.  The  velocity  mean  free  path  is 

K($)  =  vg(lq)Tv(iq)  , 

where  iv(k)  is  the  velocity  collision  time  for  the  injected  electrons.  The  velocity  relaxation 
time  reduces  to  the  momentum  relaxation  time  only  in  the  special  case  of  the  constant  effective 
mass,  which  is  far  removed  from  the  situation  here.  Figure  2  shows  the  velocity  mean  free 
path  for  several  semiconductors  as  a  function  of  energy  for  electrons  traveling  in  the  (100) 
direction.  Once  again,  the  GaosIriosAs  alloy  is  the  only  practical  material  with  a  Av  greater 
than  that  of  GaAs  (Krishnamurthy,  Sher,  and  Chen,  1988a,  1988b).  These  calculations  include 
all  scattering  mechanisms  except  the  deformation-potential-coupled  phonon  interaction.  The 
donor  concentration  is  taken  to  be  1018  cm-3.  The  T  to  L  and  X  energy  separations  in  GaAs 
arc  0.3  eV  and  0.55  eV,  respectively.  Nothing  dramatic  happens  to  the  GaAs  curve  (c)  at 
these  energies,  because  the  intcrvalley  scattering  entails  a  large  momentum  transfer,  and  the 
longitudinal  optical  phonon  and  impurity  matrix  elements  decrease  as  1/Ak  and  1/Ak2,  respec¬ 
tively.  This  is  because  these  interactions  are  long  range;  hence,  their  Fourier  components  in 
momentum  space  reach  a  peak  at  low  Ak.  While  the  deformation  potential  coupling  constant 


10 


f 

t 

► 

[■ 


E  [kV/cml 

FIGURE  1  VELOCITY-FIELD  CHARACTERISTICS  OF  GaAs-BASED  ALLOYS 


for  low  Ak  is  far  smaller  than  that  of  the  polar  optical  coupling  constant,  it  does  not  fall  for 
high  Ak  because  it  is  a  short-range  interaction.  Thus,  while  the  deformation  potential  contri¬ 
butes  little  to  intravalley  scattering,  it  is  important  in  intervalley  scattering. 

We  have  calculated  the  deformation  potential  for  intervalley  coupling  constant  for  GaAs 
from  the  underlying  band  structure  (Krishnamurthy,  Sher,  and  Chen,  1988c).  We  find  D  = 
3.25  x  108  eV/cm,  compared  to  the  generally  accepted  value  5.1  x  108  eV/cm.  Because  D2 
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FIGURE  2  VELOCITY-FIELD  CHARACTERISTICS  OF  InP-BASED  ALLOYS 


enters  into  the  scattering  range,  this  is  a  large  difference.  However,  as  discussed  in  the  intro¬ 
duction,  when  placed  in  a  full  band  structure  calculation  (where  the  density  of  states  is  higher 
than  found  in  the  effective-mass  approximation),  the  net  scattering  rates  with  the  smaller  D  and 
proper  density  of  states  and  the  larger  D  but  small  density  of  states  are  not  too  different.  We 
have  calculated  the  injection  energy  variation  of  a,  the  ratio  of  the  drain  to  source  current  for 
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\ 
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the  HET  structure.  Our  prediction  agrees  closely  with  experiment  (Krishnamurthy,  Sher,  and 
Chen,  1988c).  This  theory  has  no  adjustable  parameters.  We  conclude  that,  for  the  geometry 
used,  the  experiment  actually  reached  an  ideal  limiting  result  for  a  GaAs  HET. 

We  have  calculated  the  upper  cut-off  frequency  and  a  for  a  number  of  materials  that  are 
given  in  Table  1.  The  calculations  do  not  include  deformation-potential  scattering;  hence,  the 
reported  results  are  for  conditions  in  which  the  injection  energy  is  below  the  T-L  valley  separa¬ 
tion.  Thus,  there  is  no  intervalley  scattering.  For  the  thicker  (500-A)  devices,  InAs  and 
Gao.5lno.5As  are  better  than  GaAs.  (These  InAs  numbers  were  found  for  an  injection  energy 
below  the  bandgap  energy  to  avoid  breakdown.)  However,  for  very  thin  devices,  GaAs  is  as 
good  as  any  other  material  considered  here  and  is  the  simplest  material  to  fabricate.  Notice 
that  Gao7Al0.3As  is  substantially  poorer  than  the  others,  and  one  should  avoid  its  use  in  the 
active  regions  of  the  devices.  We  have  examined  all  the  other  semiconductors  and  none  are 
better  than  the  materials  in  Table  1  (Krishnamurthy,  Sher,  and  Chen,  1988a,  1988b). 


Table  1 

CUTOFF  FREQUENCY  vc  AND  FRACTION 
OF  INJECTED  ELECTRONS  atal 
FOR  TWO  BASE  WIDTHS 


Base  Width 
=  500  A 

Base  Width 

=  100  A 

Alloy 

RH 

Ml 

Ml 

GaAs 

3.1 

0.63 

15.6 

0.91 

Oao.5InQ.5As 

4.1 

0.82 

20.7 

0.96 

InAs 

4.3 

0.86 

21.5 

0.97 

InAso.sPo.s 

4.1 

0.76 

20.7 

0.95 

InP 

3.7 

0.64 

18.6 

0.91 

030.7^10.3  As 

1.7 

0.04 

8.3 

0.54 

13 


IV  MATERIALS  CHOICES  FOR  DRIFT-TRANSPORT-MODE  DEVICES 


A.  Electron-Temperature  Model 

The  electron  drift  velocity  v  in  many  semiconductors  increases  linearly  under  the  influence 
of  electric  field  to  reach  a  peak  value  vp  at  a  threshold  field,  ET;  any  further  increase  in  E 
results  in  reduced  v,  giving  rise  to  a  negative  differential  resistance.  High-speed  device  appli¬ 
cations  benefit  from  materials  with  maximum  vp  and  minimum  ET. 

The  v-E  curves  are  calculated  in  a  variety  of  approximations.  In  the  simplest  (two-valley) 
model,  each  valley  is  characterized  by  a  different  effective  mass.  Their  minima  are  separated 
by  an  offset  energy  AE;  these  three  parameters,  together  with  the  electron-phonon  coupling 
constant,  are  adjusted  to  fit  the  observed  v-E  curves.  Various  authors  have  evaluated  the  effect 
of  disorder  scattering  on  low  field  mobility  in  the  alloys.  The  main  difficulty  in  extending 
these  models  to  the  hot  electron  problem  is  in  applying  the  appropriate  correction  to  the  change 
in  the  effective  mass  and  scattering  potential  with  wave  vector  k.  Our  improvement  over  pre¬ 
vious  efforts  is  to  use  realistic  and  accurate  band  structures  in  evaluating  various  scattering 
rates  and  parameters  needed  to  obtain  v-E  curves  (Krishnamurthy,  Sher,  and  Chen,  1987).  In 
the  alloys,  CPA  band  structure  is  used  where  the  alloy  scattering  obtained  directly  from  the 
imaginary  part  of  the  CPA  self-energy.  In  this  formalism,  the  nonparabolic,  anisotropic  and 
intcrvalley  scattering  mediated  by  alloy  disorder  are  included  without  adjustable  parameters. 

In  this  electron-temperature  model,  the  electrons  are  assigned  a  temperature  Te,  which  is 
usually  much  higher  than  lattice  temperature  T.  This  is  because  the  average  electron  energy 
under  high  electric  field  is  much  larger  than  the  average  lattice  excitation  energy  per  particle. 
If  we  denote  the  mobility  in  the  lower  valley  as  p,  and  assume  that  the  mobility  in  the  higher 
valley  is  zero,  then: 

-l 

(4-1) 


v  =  p,E 


1  +  R  c 


AE 

kiiT. 


where  R  is  the  ratio  of  density  of  states  in  the  upper  valley  to  that  of  lower  valley.  In  the 
steady  state,  the  energy  gained  by  an  electron  from  the  field  must  be  balanced  by  the  energy 
lost  to  the  lattice: 


evE  = 


3 

2 


kB 


(Te  -  T) 

XE 


(4-2) 


where  xE  is  the  energy  relaxation  time.  Equations  (4-1)  and  (4-2)  can  be  combined  to  yield  Te 
and  v-E  characteristics.  In  the  case  of  alloys,  the  concentration  dependence  of  R,  AE  and  Pi 
are  used  to  obtain  v-E  curves.  The  results  of  these  calculations  have  been  published 
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(Krishnamurthy,  Sher,  and  Chen,  1987).  We  concluded  that  Ga0  471^  S3As  alloys  and  InP- 
based  alloys  such  as  InP075As0.25  and  InP075Sb025  are  prefera^.e  to  GaAs  for  high-speed  dev¬ 
ice  applications. 


B.  Solution  to  the  Boltzmann  Equation:  An  Eigenvalue  Method 

Although  the  full  band  structure  is  used  in  the  calculations  described  in  Section  IV-A,  the 
results  offer  only  qualitative  trends.  This  is  because  of  various  approximations  such  as  elec¬ 
tron  temperature,  displaced  Maxwellian  and  energy  relaxations,  and  the  like.  To  obtain  quanti¬ 
tatively  accurate  results,  one  must  solve  Boltzmann  transport  (BT)  equation  without  such  limit¬ 
ing  approximations.  Techniques  such  as  Monte  Carlo  combined  with  the  Rees  algorithm  are 
widely  used  to  obtain  the  distribution  function  by  solving  semiclassical  BT  equation.  The 
numerical  computation  is  tedious  and  time-consuming;  generalization  of  this  method  to  include 
full  band  structure  is  computationally  more  time-consuming.  To  circumvent  this  problem,  we 
developed  an  eigenvalue  method  to  solve  the  steady-state  BT  equation.  This  method  of  obtain¬ 
ing  solution  is  about  two  to  three  orders  of  magnitude  faster  than  the  existing  methods  (Krish- 
namurthy,  Sher,  and  Chen,  1988d). 

We  here  describe  this  method  in  brief.  The  steady  state  BT  equation  is 

-IT  ~  I  [SOcjc')  f(1c)  -  f(£',!c)  f(k*>J  =  0  .  (4-3) 

H  dkz  kv 

E  is  assumed  to  be  in  1  direction.  If  we  expand  the  distribution  function,  f(k5,  in  terms  of  a 
complete  set,  (<J),(k)},  then 


f(?)  =  I  <(>,(£)  Q 

i 


(4-4) 


Substituting  Eq.  (4-3)  in  Eq.  (4-2),  we  get  a  matrix  equation. 


MC  =  0 


(4-5) 


where  M  is  a  (n  x  n)  matrix  ,  C  is  a  column  matrix,  and  0  is  a  null  column  matrix;  n  is  the 
number  of  basis  functions  in  the  chosen  basis  set.  From  Eq.  (4-5),  we  see  that  the  eigenvector 
corresponding  to  zero  eigenvalue  of  the  M  gives  us  the  expansion  coefficients  necessary  to 
obtain  the  distribution  function,  as  given  in  Eq.  (4-4). 

If  the  chosen  basis  set  is  complete,  the  eigenvalue  zero  is  ensured  and  the  distribution 
function  can  be  accurately  obtained.  However,  if  the  basis  set  contains  very  large  number  of 
functions,  then  a  numerically  difficult  diagonalization  of  very  large  matrix  must  be  carried  out. 
Hence,  we  need  to  make  an  educated  guess  of  the  basis  set  that  is  finite  and  that  can  accurately 
describe  the  distribution  function  in  high  electric  fields.  In  all  calculations  reported  here,  we 
assumed  two  basis  function  (s  and  pz-like)  at  each  F,  L,  and  X  valleys.  The  applied  electric 
field  is  assumed  to  be  in  Z  direction.  Because  L  valleys  are  symmetrically  located  with  respect 
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to  field  direction,  the  number  of  basis  functions  in  L  valleys  can  reduced  to  just  two.  How¬ 
ever,  we  need  two  set  of  basis  function  at  X  valleys  (one  set  at  the  X  valley  in  the  2  direction 
and  the  other  set  at  X  valley  perpendicular  to  the  S  direction).  Hence,  we  have  eight  basis 
function  in  our  basis  set,  given  as  follows: 

0i  =  N,  e"“l£?  (4-6) 


02  -  ^2 
03  =  ^3 

04  =  N4 
05  =  Ns 
06  =  N6 
07  =  N7 
08  =  Ng 


deg  -“2*7 
dkz 

e-Pi(*7-tL) 

dei? 

dk2 

-Yi  «?-£,) 

-Yi  (£?-£,)  dek‘ 

dkz 


c-Y,(q*-£y) 


d£i 

dkz 


where  N,  are  the  normalization  constants.  Owing  to  the  finiteness  of  the  basis  set,  the  zero 
eigenvalue  of  M  is  not  ensured.  As  a  tradeoff  for  a  large  basis  set,  we  introduced  an  iterative 
procedure.  In  all  calculations  reported  here,  we  set  cq  =  =  a.  We  varied  a  to  obtain 

the  zero  eigenvalue.  In  most  cases,  the  zero  eigenvalue  is  obtained  in  less  than  10  iterations 
where  each  iterations  takes  approximately  2  Cray  CPU  seconds.  Once  the  eigenvectors  that 
correspond  to  zero  eigenvalue  are  obtained,  the  distribution  is  calculated  from  Eq.  (4-4).  The 
drift  velocity  can  then  be  calculated  from 

vd  =  -jjr  Z  vz(t)  f(ic)  (4-7) 


where  vz  =  dejcVdkz. 

Figures  1  and  2  show,  respectively,  calculated  v-E  curves  for  GaAs-based  alloys  and  InP- 
based  alloys.  Comparative  evaluation  is  made  in  Figure  3.  In  the  lower  panel  of  Figure  3,  Ep 
the  minimum  energy  necessary  for  pair  production  and  AE  is  plotted  for  various  alloys  and 
compounds.  Those  alloys  in  which  AE  <  Ep  cannot  be  used  efficiently  for  high-speed  device 
applications  as  the  breakdown  is  expected  before  reaching  the  peak  velocity.  The  upper  panel 
of  Figure  3  shows  vp  and  Ef 1  for  various  alloys  and  compounds.  vp  and  Ef 1  are  the  ratio  of 
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FIGURE  3  COMPARATIVE  EVALUATION  OF  VARIOUS  SEMICONDUCTOR 
ALLOYS  FOR  HIGH-SPEED  DEVICE  APPLICATIONS 


peak  velocity  and  threshold  field  in  the  material  to  that  in  GaAs.  The  materials  in  which  vp 
and  Ef 1  are  larger  than  unity,  are  preferable  to  GaAs  for  high-speed  and  power  device  applica¬ 
tions.  Our  calculations  suggest  that  Ga0.47lno  53As,  InPo.75Sbo.25>  and  InP0  75Sbo.25  are  superior 
to  GaAs. 
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V  HIGHLIGHTS  AND  CONCLUSIONS 


An  accurate  transport  calculation  for  ballistic  and  drift  mode  devices  is  reported  here.  The 
incorporation  of  this  transport  theory  into  device  modeling  codes  should  make  it  possible  to 
develop  a  computer-aided  design  model  for  these  structures.  Thus,  we  are  one  step  closer  to 
materials  engineering  to  a  purpose. 

Based  on  this  36-month  study  of  hot  electron  transport  in  semiconductors  and  alloys,  we 
conclude  the  following: 

(1)  The  v-E  curves  obtained  in  the  EMA  are  quantitatively  different  from  that  obtained 
by  using  full  band  structures  (Krishnamurthy,  Sher,  and  Chen,  1987). 

(2)  Alloy  scattering  is  important  and  affects  the  v-E  curves  substantially.  Because  this 
scattering  increases  with  density  of  states,  it  becomes  very  important  in  high  field 
transport  (Krishnamurthy,  Sher,  and  Chen,  1987). 

(3)  The  strength  of  alloy  scattering  depends  sensitively  on  the  scattering  potential  used. 
In  CPA  formalism  adapted  here,  the  scattering  potential  emerges  naturally  from  the 
underlying  band  structure  (Krishnamurthy,  Sher,  and  Chen,  1987;  Krishnamurthy, 
Sher,  and  Chen,  1988). 

(4)  Results  using  electron  temperature  two-valley  model  depend  sensitively  on  the  energy 
relaxation  time. 

(5)  Velocity  relaxation  rate  and  velocity  mean  free  path  are  the  appropriate  quantities  to 
study  the  hot-electron  transistor  characteristics  (Krishnamurthy,  Sher,  and  Chen, 
1988a;  1988b). 

(6)  The  fraction  of  electrons  that  are  ballistic  in  ballistic  devices  depend  strongly  on  elec¬ 
tron  electron  scattering.  However,  the  ratio  of  collector  current  to  emitter  current,  a, 
is  affected  by  electron-electron  scattering  only  by  10%  even  in  10,8cm-3  (Krish¬ 
namurthy,  Sher,  and  Chen,  1988c). 

(7)  The  coupling  constant  for  deformation  potential  coupled  intervalley  scattering  is  cal¬ 
culated  from  the  underlying  band  structure.  The  calculated  strength  of  the  intervalley 
scattering  explains  the  measured  a  in  GaAs  hot  electron  transistor  (HET)  very  well 
(Krishnamurthy,  Sher,  and  Chen,  1988c).  The  calculated  coupling  constant  is  much 
smaller  than  the  ones  currently  used  in  the  literature.  Because  the  measured  a  agrees 
with  our  calculated  one,  it  demonstrates  that  a  near-ideal  HET  was  made  and  no 
further  increase  in  its  function  should  be  expected. 

(8)  Owing  to  the  changing  curvature  of  the  band  structure,  the  obtained  velocity  scatter¬ 
ing  rates  for  impurities  and  phonons  are  larger  than  that  obtained  in  EMA. 

(9)  For  ballistic  transport  devices  with  thicker  width  (approximately  500  A),  InAs  and 
Ga^In^As  materials  arc  expected  to  be  more  efficient  than  GaAs  (Krishnamurthy, 
Sher,  and  Chen,  1988a;  1988b). 


(10)  For  ballistic  devices  with  very  thin  base  width  (approximately  100  A),  no  material 
offers  greater  advantages  than  GaAs  (Krishnamurthy,  Sher,  and  Chen,  1988a;  1988b). 

(11)  We  had  developed  a  new  and  faster  numerical  method  to  solve  the  steady-state 
Boltzmann  transport  (BT)  equation.  This  iterative  solution  to  BT  equation  can  be 
obtained  in  typically  2  cpu  seconds  on  the  CRAY-2  (Krishnamurthy,  Sher,  and  Chen, 
1988d). 

(12)  Using  only  eight  basis  functions,  we  have  obtained  v-E  curves  up  to  15  kV/cm  of 
electric  field.  The  calculated  peak  velocity  and  threshold  field  are  in  excellent  agree¬ 
ment  with  experiments.  It  should  be  pointed  out  that  these  results  are  obtained  with 
absolutely  no  adjustable  parameter. 

(13)  Our  calculations  suggest  that  Gao.47Ino.s3 As,  InP0  75As0.25 As  and  InP0  75Sbo.2j As  are 
better  than  GaAs  for  high-speed  device  applications.  Gao.47Ino.33As  is  the  most 
efficient  of  all  the  candidates  considered  here  for  these  applications  (Krishnamurthy, 
Sher,  and  Chen,  1988d). 

(14)  To  obtain  solution  to  BT  equation  at  higher  electric  fields,  the  number  of  basis  func¬ 
tions  must  be  increased  to  describe  the  distribution  function  adequately.  This 
improvement  should  be  carried  out  to  obtain  peak-to-valley  ratio  and  saturation  drift 
velocity  in  these  alloys. 

(15)  This  eigenvalue  method  should  be  slightly  modified  for  the  cases  of  spatially  varying 
electric  field. 
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